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We confirm by explicit computation the conjectured all-orders iteration of planar maximally su- 
persymmetric A/" = 4 Yang-Mills theory in the nontrivial case of five-point two-loop amplitudes. We 
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using a Mellin-Barnes representation and the automated package of ref. [1]. This confirmation of 
the iteration relation provides further evidence suggesting that A/" = 4 gauge theory is solvable. 
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In his seminal work dating to the infancy of asymp- 
totic freedom, 't Hooft [2] gave hope of solving quan- 
tum chromodynamics (QCD) in the so-called planar 
limit, when the number of colors is taken to be large. 
While this hope for ordinary QCD has not yet been real- 
ized, the Maldacena conjecture [3] has brought it closer 
for four-dimensional maximally-supersymmetric Yang- 
Mills theory (MSYM), by proposing a duality relating 
it at strong coupling to type IIB string theory in five- 
dimensional anti-de Sitter (AdS) space at weak cou- 
pling. Heuristically, this suggests that the leading-color 
terms of the perturbative series should be resummable, 
and along with possible non-perturbative contributions 
should yield relatively simple results matching those of 
weakly-coupled gravity. 

While the Maldacena conjecture does not address di- 
rectly the scattering amplitudes of on-shell (massless) 
quanta, previous work by Anastasiou, Dixon, and two 
of the authors [4] shows that the basic intuition holds. 
That paper presented a conjecture for an all-orders iter- 
ative structure in dimensionally-regulated scattering am- 
plitudes of MSYM. Dixon and two of the authors [5] 
fleshed out this structure for maximally helicity-violating 
(MHV) amplitudes. Witten's proposal [6] of a weak- 
weak duality between MSYM scattering amplitudes and 
a twistor string theory provides further indications of new 
structures underlying the simplicity of both MSYM and 
string theory in AdS space at strong world sheet coupling. 

Ref. [4] verified the iteration conjecture explicitly for 
the two-loop four-point function (a second verification 
was given in ref. [7]), and ref. [5] did so for the three-loop 
four-point amplitude. Furthermore, the computation of 
the two- loop splitting amplitude in ref. [4] , its own itera- 
tion relation, and consideration of limits as momenta be- 
come collinear shows that were the conjecture to hold for 
the five-point two-loop amplitude, it would almost cer- 
tainly hold for all MHV two- loop amplitudes. The step 



from four-point to five-point amplitudes is non-trivial, 
because at five points, functions that are not detectable 
in the collinear limits appear [8]. 

An important step in closing this gap has recently been 
taken by Cachazo, Spradlin and Volovich [9]. They con- 
firmed the conjecture for the terms in the two-loop five- 
point amplitude even under parity, using an earlier guess 
for the integrand [10]. In this Letter, we will complete 
the task. We compute the integrand using the unitarity 
method [8, 11, 12], confirming the form of ref. [10] for 
the parity-even terms, and providing the correct form for 
the parity-odd ones. We then integrate numerically at 
random kinematic points, using the MB integration pack- 
age [1], to show that the conjecture holds for both parity- 
even and -odd terms. We also remark that the 'addi- 
tional iterative structure' of ref. [9] follows from the one 
of ref. [4] by setting odd parity terms to zero on both 
sides of the iteration formula. 

The unitarity method [8, 11, 13] has proven powerful 
for computing scattering amplitudes of phenomenological 
and theoretical interest out of reach using conventional 
Feynman diagrammatic methods. Improvements [14] 
have followed from use of complex momenta [6] . 

Perturbative amplitudes in four-dimensional massless 
gauge theories contain infrared singularities. These are 
well understood [15] in MSYM and are a subset of the 
ones appearing in QCD. As in perturbative QCD, the S- 
matrix under discussion here is not the textbook one for 
the 'true' asymptotic states of the four-dimensional the- 
ory, but rather for states with definite parton number. As 
in QCD, a summation over degenerate states would be re- 
quired to obtain finite results for scattering [16]. We reg- 
ulate these divergences in a supersymmetry-preserving 
fashion using the four-dimensional helicity (FDH) [17] 
variant of dimensional regularization, with D ^ A — 2e. 
(This scheme is a close relative of Siegel's dimensional 
reduction [18].) 
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FIG. 1: The three- and two-particle cuts of the five-point 
amphtude. 



We write the leading-color contributions to the L-loop 
SU{Nc) gauge-theory n-point amplitudes as, 
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where 7 is Euler's constant, the sum is over non-cyclic 
permutations of the external legs. We have suppressed 
the momenta and helicities ki and A^, leaving only the 
index i as a label. This decomposition holds for all parti- 
cles in the gauge super-multiplct as all are in the adjoint 
representation. We will find it convenient to scale out 
the tree amplitude, defining MA (e) = An /An ■ 

At two loops the iteration conjecture expresses n-point 
amplitudes entirely in terms of one-loop amplitudes and 
a set of constants [13]. For MHV amplitudes up to 0(e'^), 

Mi'Ke) ^ liMi'He))' + f^'\e) M(^\2e) + C^'^ , (2) 

where /(2)(e) = -(^^+(^e+C4e^+- ■ •), and C(2) = -(1/2. 
Ref. [5] provides analogous equations for higher-loop 
MHV amplitudes. Subtracting out the known infrared 
divergences [15] provides an all-loop form for the finite 
remainder, expressed in terms of the one-loop finite re- 



mainder and two constants, one of which is an anoma- 
lous dimension. A conjecture for the required anomalous 
dimension was very recently presented [19], based on a 
proposed all- loop Bethe Ansatz [20]. It is rather inter- 
esting that this anomalous dimension corresponds to one 
of the terms appearing in the QCD one [21]. 

To check whether the iteration relation holds in the 
critical five-point case, we have evaluated a set of cuts suf- 
ficient to determine the five-gluon integrand completely. 
These include the thrce-particlc cuts depicted in fig. la 
as well as the contributions to the two particle cuts from 
fig. lb. The three-particle cuts on their own determine 
all integral functions, except for those which are simple 
products of one-loop integrals. The two-particle cuts rule 
out the latter (double cuts suffice). 

The use of a dimensional regulator involves an analytic 
continuation of the loop momenta to D dimensions. At 
one loop, the discrepancy between treating loop momenta 
in four or D dimensions does not modify the amplitudes 
of a supcrsymmetric gauge theory through 0(e'^). No 
such proof exists for higher loops. Thus, to ensure that no 
contributions are dropped, we compute the unitarity cuts 
in D dimensions [22]. This does complicate the analysis, 
because standard helicity states can no longer be used 
as the intermediate states. We can avoid some of the 
additional complexity by considering instead the D = 
10, A/" = 1 super- Yang-Mills theory. When compactified 
on a torus to Z? = 4 — 2e dimensions this is equivalent to 
dimensionally-regulated MSYM in the FDH scheme. 

After reducing all tensor integrals we obtain an ex- 
pression for the amplitude in terms of the integrals 
shown in fig. 2. The color-ordered amplitude with four- 
dimensional external momenta is given by a sum over the 
cyclic permutations of those momenta, 
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Here, s^ = {h + kjY , R = £1234512523534545551/61 



234, 



Sabc = 512S51 + aSi2523 + ^523534 - S51S45 + CS34S45 , 

£1234 = Aie^.^p^k'^kl^kf^kl = tr[75^i/^2fe^4] , (4) 

and G1234 — det(sy), (i,j ~ 1,...,4). (In 5, a,b,c ~ 
±1.) The terms lacking a factor of £1234 are even under 
parity, while those with such a factor are odd. The even 
terms match the guess originally given in ref. [10], but 
the odd terms differ (the odd terms in ref. [10] do match 



the four-dimensional double two-particle cuts). 



Because of the l/e^ infrared singularity in one-loop am- 
plitudes, and because these appear squared in the itera- 
tion relation, we need expressions valid through 0{e^). A 
representation of the one-loop five-point amplitude, ex- 
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FIG. 2: The two-loop integrals appearing in the five-point 
amplitude, with all external momenta flowing outwards. The 
normalization is as given in eq. (8), and the numerical labels 
on the internal propagators in (c) specify the arbitrary powers 
fli. The prefactor in (c) is understood to be inserted in the 
numerator with power — ag; in eq. (3), — ag = 1. 
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FIG. 3: The one- loop integrals required to all orders in e 
for the one-loop five-point amplitude. The normalization is 
as given in eq. (7), and the numerical labels on the internal 
propagators in (b) specify the arbitrary powers a^. 



tending ref. [8] to all orders in e, may be found in ref. [23], 
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in terms of the integrals of fig. 3. As indicated by the su- 
perscript, the second integral (fig. 3b) is to be evaluated 
in 6 — 2e dimensions. In D = 6 it is completely finite, 
but because it appears multiplied by an infrared-singular 
integral in eq. (2) we need its value through 0{e). 

To obtain Laurent expansions in e for our integrals 
we use the Mellin-Barnes (MB) technique, successfully 
applied in numerous calculations (see, e.g., refs. [5, 24, 
25, 26, 27] and chap. 4 of ref. [28]). It relies on the 
identity 
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where —Re A < /3 < 0. This basically replaces a sum over 
terms raised to some power with a product of factors. 

The box function in fig. 3a was given to all orders in e in 
terms of a hypergeometric function in ref. [29] . Here we 
need its value through 0{£^). Evaluating the pentagon in 
fig. 3b with arbitrary powers of propagators also allows a 
parallel evaluation of this integral to the required order. 

The derivation of a fourfold MB representation for 
the one-loop pentagon diagram is straightforward, after 
Feynman parameterizing. 
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where Kij ^ h + kj , 02345 = a2 + a^ + a4 + a5,A = 
Yl o,ij ^124 = zi + Z2 + Z4, etc. We have allowed for ar- 
bitrary powers of propagators so that we can obtain all 
one-loop integrals. Taking 05 — ^ 0, with other a; = 1, 
gives the box integral /K in fig. 3a. Setting all a, == 1 
and shifting all terms except the e'^'^ prefactor by e — > e— 1 



J 



yields the D = 6 — 2e pentagon, I„! '^, corresponding 
to fig. 3b. 

The various two-loop pentabox integrals have a sev- 
enfold MB representation obtained by inserting a three- 
fold MB representation for a two-mass double box into 
eq. (7), 
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The limit ag ^ or 07 ^ with ag = and the 
other a,i = 1 yields the double box with one massive leg 
(fig. 2a and b) in agreement with ref. [27, 30]. Moreover, 
P(2)(i^ . . . , 1^ _1)^ p(2)(i^ . . . , 1,0) and P(2)(i^ ... 1,0,0) 
yield the integrals in fig. 2c, fig. 2e and fig. 2d, respec- 
tively. 

An essential step in the use of the MB technique is 
the resolution of singularities in e or zeros that appear as 
tti — > 0. There are two strategies for doing this [24, 25]. 
Quite recently, the second strategy was formulated algo- 
rithmically [1, 31] and implemented in the MB package [1]. 
It produces code that allows the integrals to be evaluated 
numerically to reasonably high accuracy. 

The even terms in eq. (3) were recently evaluated in 
ref. [9] using the MB package [1] along with the guess of 
ref. [10]. Those authors observed that an iterative struc- 
ture holds for the parity-even terms alone. We may ob- 
serve that this structure is not independent of the com- 
plete iteration formula eq. (2): use the results for the one- 
loop five-point amplitude in eq. (5), set the odd terms to 
zero and use the fact that the one-loop MHV amplitudes 
have even parity through 0{e^). At higher loops, we do 
not expect a clean separation between 'even' and 'odd' 
terms, as non- vanishing terms of the form £^234 ^il^ arise. 
These are even under parity. 

We have evaluated all the two-loop integrals in fig. 2 
through 0{e^^) and the one-loop integrals in fig. 3 through 
©(e^) using the representations in cqs. (7) and (8). We 
have checked to a numerical accuracy of five significant 
digits at three independent kinematic points that the it- 
eration formula (2) is indeed correct for the complete am- 
plitude. This is a crucial check on the conjecture because 
the parity-odd terms in the five-point amplitude are pre- 
cisely the ones which are not constrained by coUinear 
factorization onto four-point amplitudes. 

The calculation presented here makes a nontrivial ad- 
dition to the existing body of evidence for the iteration 
conjecture [4, 5]. The conjecture itself gives us good 
reasons to believe that MSYM is solvable. Within the 
context of the planar perturbative expansion, this would 
imply the resummability of the series. Parallel develop- 
ments in uncovering the integrable structure of the theory 
(see e.g. refs. [20, 32]) also lend credence to this belief. 
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